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ABSTRACT 

The wave functions in the diquark- quark model for baryons are modified to take 
into account the effect of the different masses of the quarks in diquark formation. A 
spatial wave function is introduced multiplying each flavour-spin configuration in the 
quark-diquark model. Futhermore a numerical coefficient is introduced which weighes 
the contribution of each configuration. 
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1. Introduction 



The concept of diquark was introduced, in the beginnings of the quark model to 
describe baryonic properties. Diquarks may be conceived as substructures or clusters 
of pairs of quarks inside baryon wave function [1,2]. Diquark-quark models have been 
intensively developed in many works on baryon structure and mass spectroscopy, baryon 
decay and high energy baryon scattering [3]. Some diquark-quark wave functions were 
constructed to describe the structure of baryons. Those wave functions have the disad- 
vantage of not being completely antisymmetric under exchange of any two quarks. In 
fact, they are not antisymmetric under the exchange of a quark of the diquark and the 
third quark. This dificulty is compensated by the simplicity of the model. However, the 
absence of antisimmetry is consistent with the diquark interpretation as an quasielemen- 
tary object. Futhermore in the baryon wave function construction [4] mentioned above, 
it was not taken into account the effect of the difference in the masses of quarks. 

The interactions among the quarks inside the baryons are independent of the 
flavours, but the diquark dimension depends on the masses of the constituent quarks. 
The higher the quark masses the smaller the dimension of the diquark. 

Given a baryon with three quarks of different masses, there are three different 
configurations of diquark-quark that can contribute to its wave function. The diquarks 
in each configuration are not equally probable. Another serious distinction in the three 
possible configurations of the diquark-quark model is that their masses are different. 
This result constitutes an ambiguity in the diquark-quark model for baryon. 

In this paper we give a prescription to construct the baryon wave function, in 
which the contribution of each diquark-quark configuration is adequately taken into 
account. The expectation value of the Hamiltonian calculated with the wave function 
obtained by this prescription gives the baryon mass. The mass of the baryon is, then, 
the weighed mean value of the masses of the diquark-quark configurations. This mean 
value eliminates the ambiguity in the computation of the baryon mass. The prescription 
is applied to the construction of the wave function of the lowest energy baryonic states 
in S wave with spin 1/2. In these wave functions there exists two spin states, S=0 e S=l, 
for each S wave diquark. Then each flavour configuration has a scalar and a vectorial 
diquark. 

2. Kinematics and Diquark approximation 

The coordinates used in the description of the motion of the three quarks in the 
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center of mass system of the baryon are 



Pi = Tj - f k (1) 

and 

- _ rrijfj + m k f k 

K = ri LJ — (2) 

rrij + m k 

The conjugate momenta of these coordinates are 

Pi = HjkPi (3) 

and 

qi = fiiXi, (4) 

where 

Pjk = (5) 

rrij + m k 

mi(mj+m k ) 

= (6) 

mi+m 2 + m 3 
The kinetic energy of the three quarks is 

T = + ^- (7) 

In the diquark model, the potential energy can be written, as 

V*V jk (pi) + Vi dk (\i). (8) 

This is an approximate form of the potential, consistent with the assumption that 
diquarks are quasielementary objects inside the baryon. 

The Hamiltonian of the baryon, in the diquark model, decouples into two Hamil- 
tonians, 

H » + H ijjk (9) 
The internal motion of the diquark(jk) is described by 

Hjk=7T- + V MP^ ( 10 ) 
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and the Hamiltonian of the relative motion of the diquark(jk) and the third quark(i) is 

H i,jk = + v i,jk{K) (11) 

The wave function of the three quarks for the i(jk) quark-diquark configuration 
factorizes as 



i/>i(pi, X i) = <t>jk(pi)Xi(h) (12) 
Assuming that 4>jk and Xi are respectively eigenfunctions of 

Hjk4>jk = E jk (f)j k (13) 

and 

HijkXi = EijkXi (14) 
The energy of the configuration i(jk) is given by 

Ei = Ej k + Eijk (15) 
3. Potentials in diquark model 

The potentials Vjk(pi) and Vijk(K) have the general form 

V(f) = V{r) + U apin {f) (16) 
The spherically symmetric part of the potential can be written as 

V(r) = C + V v (r) + V s (r) (17) 

V v (r) = V coul (r) + (1 - f)V conf (r), (18) 

where 

V s (r) = fV conf (r) (19) 



Vooul = —Fq— (20) 
r 

and 



V conf = Kr* 



(21) 



The parameter f controls the proportion of scalar and vectorial parts under Lorentz 
transformation of the confining potential [5,6]. 
For each pair(ij) of particles, Fq is given by 

F G = (F i .F j ) = F? j -F?-F] (22) 

where F^ is the expectation value of the quadratic Casimir operator for the colour group 
SU(3) in the (i) multiplet. 

The parameter C depends on the masses of the interacting particles (ij) according 

to 



C(rrii, rrij) = ao + a\x + a^x , (23) 

where 



x = ln(mfmj + rajm^)). (24) 

The running coupling constant of chromodynamics, a s , and the parameter K of the 
confining potential are assumed to be constant for a given colour system and independent 
of the masses of the interacting particles. 

The potential defined by eqs. (14) to (21) was applied to the study of the spec- 
troscopy of the mesons (qq) [6], glueballs [7], and hybrids (qqg) [8]. The masses of the 
quarks and the parameters of the potencial used in this work are shown in Table [1]. 

The spin-dependent terms of the potential can be obtained from the Breit-Fermi 
Hamiltonian, [9]. 

In this paper, only S-wave states are computed. Then the only contributing spin 
dependent term is 

V ss = (-^— ) V 2 K(^ • S 2 ). (25) 

The lowest energy states can be obtained by a variational approach, choosing gaus- 
sian wave functions 

9 3 / 4 1 r 2 

0(r) = (-) exp[-(-)}. (26) 

n TV 'o 

The energy and radius po of the diquark(jk) are given by minimizing 
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Ejk < < 4>jk I Hj k | 4>j k > ■ (27) 

Similarly, the energy and radius A D of the configuration i(jk) are obtained by min- 
imizing 

Ei,jk < <Xt\ Hijk | Xi > ■ ( 28 ) 
In tables (1) to (4), we show the results of these computations. 

4. The baryon wavefunctions in quark-diquark model 

Diquarks are assumed to be antitriplet colour states, which are formed by two 
colour triplets, the quarks interacting attractively [10]. 

The flavour-spin wave functions of the diquarks are denoted by: 

V +1 {q'q) = ±={qq' + q'q)®{]]) (29a) 
V +1 (q'q) = ±=(qj + q'q) ® ^=(U + IT) (296) 
V (q'q) = -L(qq' + q'q)®ai) (29c) 
S(q'q) = -±=(qq' - q'q) ® -J=(|T " T!) (29d) 

where q,q'=(u,d,s,c,b). 

The baryon wave functions in the quark-diquark model [4] can be modifyed to take 
into account the effect of the difference in quark masses. This is done by introducing 
in the baryon wave functions [11] a factor multiplying each quark-diquark configuration 
QaiQbQc): which distinguishes among the contributions of each configuration . These 
factors are given by 

Pbc 



where (be) is a scalar diquark, and 



& = > ( 31 ) 

Pbc 
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where (be) is a vectorial diquark; p s hc and p v hc are the r.m.s radii of the scalar and 
vectorial(bc) diquarks, respectively. The r.m.s radii of the configuration a(bc) are A„ 
and A„ for scalar and vectorial diquark (be) , respectively. The parameters v a are fixed 
by imposing orthogonality on the baryon wave functions. The spatial wave functions of 
a quark-diquark flavour configuration a(bc) will be denoted by 4> s a or (jf a , according as 
the diquark(bc) spin is or 1. 

The wave functions for baryons in the flavour SU(3) octet, (uds), in the quark- 
diquark model, are given by 



p, ± >= ±F p { (3 u [V (ud)u± - V2V±(ud)u T }(f> v u + p d [-V2V (uu)d± + 2V±(uu)d T }(f> v d + 



±3>a u S (ud)u±4> 



(32a) 



a, ± >= ±F n \ p d [V (du)d± - ^V ± {du)d T ]^ v d + p u [-V2V {dd)u± + 2V±(m/)d T ]C+ 

±3a d S(du)d±(f> d ^ (326) 



-V2V±(us)u T }<Pl+P s [-V2V (uu)s ± +2V ± (uu)s T }<j) v s + 



±3a u S(us)u±(j) 



(32c) 



± >= ±F E - \ P d [V {ds)d±- V2V±{ds)d T ](t> v d +p s [-V2V (dd)s± + 2V±{dd)s T ](t> v s + 

±3a d S(ds)d±(p d \ (32d) 



S°,± >= ±F sO ^0 u [V o (ds)u± - V2V ± (ds)u T ](f>l + Pd[V (us)d± - V2V±(us)d T ](j) d - 
-2(3 s [V {ud)s± - V2V±(ud)s T ] +<P V S ± 3[a d S(su)d±<j) d + a u S(ds)u±<j) s u ] \ (32e) 
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A , ± >= ±F A0 ^p d [-V (us)d± - \/2V±{us)d T \(j) v d + P u [V (ds)u± - V2V ± (ds)u T }cf> v u + 
T[a u S{ds)u±(t> s u + a d S(us)d ± cf> d - 2a s S(du)s ± c/ ) s s ] j (32/) 

3°, ± >= ±F E o [p a [Vo(us)s± - V2V±(us)stW s +[3 u [-V2Vo(ss)u± + 2V±(ss)u t W u + 

±3a s S(us)s ± r s } (32(?) 



S",±>= ±F s -{&[Vb(<fe)s±-V^±^ 

±3a 8 S(ds)s±(j) s s \ (S2h) 



The wave functions of all baryons denoted by (A,B,C) = (u,d,s,c,b) can be written in any 
one of the general forms: 



AAB,±> 1 = ±F 1 ^p A [V (AB)A± - V2V±(AB)A-_ 
+(3 B {-V2V (AA)B± + 2V±{AA)B T ](f> v B + 

±3a A S(AB)A± ( f> s A ^ (33a) 

ABC, ± > 2 = ±F 2 ^p A [V (BC)A± - V2V±{BC)A T ]<t> v A 

+Pb[V (AC)B± - V2V±(AC)B T ]cP B + 
-2p c [V (AB)C± - V2V ± (AB)C T }<t> v c 

±3[a A S(BC)A±(j) A + a B S(AC)B±)](j) B j (336) 

ABC,± > 3 = ±F 3 ^p A [V (BC)A± - V2V±{BC)A T ]<t> v A 
+Pb[-V (AC)B ± - y/2V±(AC)B T ](f> v B 



T[a A S(CB)A±(f) s A + a B S(CA)B±(j) s B - 2a c S(AB)C±(f) a c } j. (33c) 



The normalization constants are 

32 \ -1/2 



Fi = (9c£ + 3p 2 A + 6p%) 1 (34a) 

/ \-l/2 

F 2 = \9(a 2 A + a%) + 4/£)j (346) 

F 3 = [a 2 A + a 2 B + Aa 2 c + S(p 2 A + (3 2 B ) \ (34c) 
The orthogonality conditions, which determine the values of the parameters v a , are 

2 < ABC, ± | ABC, ± > 3 = 0, (35) 



which gives 



W 2 A -pl)-(a 2 A -a 2 B ) = (36) 



The values of the parameters v a satisfying these conditions are given in Table [2]. 
The values of v a for states like | AAB, ± >i are undetermined by these conditions, and 
are assumed to be v a = 0.6, a mean value in the interval [0.2,1.0] that contains the 
values of v a shown in Table [2]. Choosing v a = 0.6, we have a general value that givesthe 
lowest values for the masses of states like | AAB >i, which is consistent with a linear 
variation function approach for baryon wave functions. 

Using the above defined notations (32), the wave functions of the spin 1/2 charmed 
baryons pertaining to 20-plet of the flavour SU(4) can be written as in Table [3]. Simi- 
larly, the wave functions of the spin 1/2 baryons with bottom flavour are shown in Table 
[4]. The values of the factors a and (3 defined in eqs (30,31), multiplying each flavour 
configuration a(bc), are listed in Tables [5] and [6]. 

The differences in the values of the a a ,ab,cx c and [3 a ,[3b,flc for a baryon with 
flavours(a,b,c) show that the formation of the diquarks (bc),(ac) and (ab) are not equally 
probable inside the baryon. The most probable configuration is the one in which the 
lighest quark is outside the diquark. 

The mass of the baryon with quarks(A,B,C) is given by 

M( A ,b,c) =m a + m b + m c + < ABC | H \ ABC > . (37) 
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The expectation values of the Hamiltonian for each one of the states (33) are 



E(AAB) 1 = Fl (30a E {AB)A + §13b E IaA)B + 9 a A E (AB)A) ( 38a ) 
E (ABC) 2 = F 2 (^flA E (BC)A+3f3B E (AC)B + 12 Pc E (AB)C + 9( a A E 

(386) 



E 



(ABC)s = E i (^(^A E IbC)A + 3Pb E (AC)B + a A E (BC)A + a B E (BC)A + 4a C'- E '(AB)C^ 

(38c) 

The masses of the spin 1/2 S-wave baryons obtained by our method are shown in 
Tables [7] and [8], and compared with the existing experimental values. 
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Table [l]-values of parameters 

m u = md = 0.38GeV", m s = 0.5GeV, m c = l.hGeV, m b = A.hGeV 
a Q = 0.0010, ai = 0.146, a 2 = -1.412 
/ = 0.5, a s = 0.187, K = 0.767 
Fq = —2/3(diquark),FG = —4/3(meson) 



Table [2] -Values of the parameters u a , for the 
quark-diquark configurations a(bc). Where (q=u,d). 
Bar. v q v s v c v h 

qsc 0.6957 0.9824 0.5 

qsb 0.3934 0.5 0.4443 

qcb 0.7811 0.5 0.9823 

scb - 0.3472 0.5 0.2493 
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Table[3]-Wave functions, in the quark-diquark model, for 
baryons with charm flavour, in the spin 1/2 S wave. 



Wave Functions with charm quark 





UUC, ± >i : 


; |s+ ±> = 


-udc, ± > 2 


|s°,±> = 


ddc, ± >i ; 


|A c h ,±> = 


udc, ± > 3 


S+±> = 


use, ± > 2 ; 


|S+*,±> = 


W.SC, ± > 3 


|s°,±> = 


dsc, ± >2 ; 


|S°*,±> = 


dsc, ± >3 


|n°,±> = 


ssc, ± >i ; 


\n+,±> = 


CCS, ± >i 


s+ + ,±> = 


ecu, ± >i : 


■ 1 H+ ± > = 


ccrf, ± >i 



Table [4] -Wave functions, in the quark-diquark model, for 
baryons with bottom flavour, in the spin 1/2 S wave. 



Wave Functions with bottom quark 



ls+±> = 


-u/u6, ± >i ; 


|S 6 -,±> = 


wd6, ± >2 


ls 6 -,±> = 


ddb, ± >i ; 


|Ag,±> = 


wd6, ± >3 


|sg,±> = 


ws6, ± > 2 ; | 


S°*,±> = 


usb, ± >3 


|5 6 -,±> = 


<is6, ± >2 ; | 


H"*,±> = 


ds6, ± >3 


|O b ",±> = 


ss6, ± >i ; 


|n»,±> = 


s66, ± >i 


|3&,±> = 


66u, ± >i ; 


l s ^±> = 


d66, ± >i 


l^,±> = 


ucb, ± > 2 ; 


"cb ) 31 ^ — 


-uc6, ± > 3 


1 =° + > - 


<ic6, ± >2 ; 


1 ^0* _i_ - _ 
1 "cb > 31 > — 


deb, ± > 3 


=° + > - 


sc6, ± > 2 ; 


1 "sbb> 31 ^ — 


| sc6, ± > 3 


n»c»±> = 


66c, ± >i ; 


l^±> = 


cc6, ± >i 
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Table [5]-Values of the factors a a for configuration 



a(bc) with scalar diquark (be). Where (q=u,d). 



diquark 


Otq 




a s 


a c 


OLh 


qq 


0.8443 


0. 


.7818 


0.6891 


0.6468 


qs 


0.8508 


0. 


.7860 


0.7263 


0.7127 


ss 


0.8664 






0.6855 


0.6301 


qc 


0.8664 


0, 


.6779 


0.6394 


0.3824 


sc 


0.8655 


0. 


.8571 


0.6713 


0.7979 


cc 


1.0358 


0. 


.9825 




0.6438 


qb 


0.8627 


0. 


.8153 


0.6550 


0.4809 


sb 


0.9290 


0. 


.8655 


0.6909 


0.5119 


cb 


1.1864 


1. 


.0369 


0.8412 


0.6394 


bb 


1.3091 


1. 


,2346 


0.9754 





Table [6]-Values of the factors (3 a for configuration 



a(bc) with scalar diquark (be). Where (q=u,d). 



diquark 


P q 








Pb 


qq 


0.4411 





.6776 


0.5927 


0.5601 


qs 


0.7719 





.6966 


0.6542 


0.6531 


ss 


0.7964 






0.6152 


0.5702 


qc 


0.8117 





.5980 


0.5853 


0.3375 


sc 


0.8047 





.8114 


0.6224 


0.7725 


cc 


1.001 





.9477 




0.6111 


qb 


0.8217 





.7758 


0.5943 


0.3937 


sb 


0.9011 





.8249 


0.6448 


0.4457 


cb 


1.1454 


1 


.0199 


0.8145 


0.6061 


bb 


1.2861 


1 


.2122 


0.9530 
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Table [7]-Masses (in GeV) of S wave and spin 1/2 baryons 



with experimental datas [12]. 



Flavour 


Baryon 


Our Results 


Experimental 


uud 


P 


1.0402 


0.9383 


udd 


n 


1.0402 


0.9396 


uus 


E+ 


1.2064 


1.1894 


uds 


E° 


1.2064 


1.1974 


uds 


A 


1.1823 


1.1156 


dds 


E- 


1.2064 


1.1926 


uss 


S° 


1.3381 


1.3149 


dss 




1.3381 


1.3213 


uuc 


E++ 


2.5193 


2.453 


udc 


s+ 


2.5193 


2.453± 0.003 


udc 


A c + 


2.4296 


2.285 


ddc 




2.5193 


2.453 


udb 




5.8639 


5.641± 0.050 
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Table [8]-Masses (in GeV) of S wave and spin 



1/2 baryons without experimental datas. 



Flavour Baryon Our Results 



use 


S + £ + * 

c ' l—l c 


2.6534,2.5843 


dsc 


= 77O* 

C ' 1 'c 


2,6534,2.5843 


ssc 


c 


2.7536 


ucc 


1 'cc 


3.7580 


dec 


s+ 

1 'cc 


3.7580 


ssc 


cc 


3.8607 


uub 





5.9862 


ddb 


s b~ 




5.9862 


usb 


770 770* 


6.1154,6.0214 


dsb 


"b ' "b 


6.1154,6.0214 


ssb 


o- 


6.1948 


ucb 


77+ 77+* 

"cb' "cb 


7.0949,7.1241 


deb 


770 770* 
1—l cb' 1 cb 


7.0949,7.1241 


scb 


770 770* 
"scb' "scb 


7.1719,7.1997 


ccb 


^c + cb 


8.2060 


ubb 


770 
"bb 


10.3836 


dbb 


"bb 


10.3836 


sbb 


^b"b 


10.4961 


ebb 


^bc 


10.4961 
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